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Both electric and magnetic screening masses in a nonperturbative gluonic background are investi- 
gated using operator product expansion. The magnetic screening mass is found to agree with lattice 
results whereas the electric screening mass is somewhat smaller than the one found on the lattice. 



PACS numbers: 11.10.Wx,12.38.Lg,12.38.Mh,25.75.Nq 

The long-range properties of a thermal system are 
characterized by screening masses or the inverse of equal- 
time correlation lengths. For a given theory the screening 
masses determine the infrared sensitivity of various ther- 
modynamic quantities as well as the spectral properties 
of the system [1 . In QED the screening associated with 
the electric fields exhibits a non-vanishing Debye mass 
whereas that with the magnetic fields does not show up 
due to gauge invariance. This indicates that the longest 
length scale in hot QED plasma is dominated by mag- 
netic fields. On the other hand, in QCD the scenario 
is far more complicated than in QED due to the gauge 
dependence of the chromo-electric and magnetic fields, 
which leads to subtleties in the calculations. 

The structure of QCD, at least near to phase tran- 
sition, seems more complex than one usually expects. 
Perturbative predictions are upset by the presence of 
strong non-perturbative effects [2 . The nonperturba- 
tive determination of screening masses has been per- 
formed in lattice QCD (LQCD ) 0-[5] with an appro- 
priate gauge fixing. The data are consistent with an over 
all exponential behaviour for the electric screening func- 
tion in all gauges whereas that for the magnetic sector 
involves a nontrivial behaviour [3, 6 . In perturbation 
theory the electric screening mass to the lowest order 
is obtained as m^) ^ gT (the strong coupling constant 
is as = ^^/47r and T is the temperature), which falls 
short of the nonperturbative description. On the other 
hand, magnetic screening cannot be addressed in per- 
turbation theory but one expects the magnetic mass to 
be generated (m^ ^ gmn ~ g'^T) nonperturbatively in 
the static sector OIH]. Nonetheless, perturbative meth- 
ods could only be reliable for temperatures far above the 
critical temperature. Moreover, the perturbative power- 
counting hierarchy of scales tud > Trim ^ 9^T^ is doubt- 
ful close to the critical temperature. Effective models 
using Polyakov Loop correlation [9 , dimensionally re- 
duced QCD flO], A/" = 4 supersymmetric Yang-Mills the- 
ory [11 and AdS/QCD [12 were employed to analyze 
the screening masses from gauge invariant correlators. 



We note that the contributions from the nonperturba- 
tive magnetic sector reveal a strong dependence on these 
correlation functions but provide very useful information 
after all. 

It is widely accepted that the nonperturbative dynam- 
ics of QCD is signaled by the emergence of power cor- 
rections in physical observables. These nonperturbative 
corrections are introduced via non-vanishing vacuum ex- 
pectation values of local quark and gluonic operators such 
as the quark condensate (V^V^) and the gluon condensate 
which are also measured in LQCD [13]. This 
approach of the operator product expansion (OPE) has 
met noticeable success in QCD sum rule calculations at 
zero temperature [14] and the calculation of the A^-point 
functions in QCD at zero temperature [15]. Unlike QCD 
sum rules, the condensates do not appear in a gauge in- 
variant combination in QCD Green's Function [16 . In 
addition, there is an explicit dependence on the gauge fix- 
ing parameter (^) in the Wilson coefficients. The OPE 
has also provided some insight on the nonperturbative 
features of QCD at finite temperature. In-medium non- 
perturbative chiral quark propagator [17^, quark-photon 
vertex [18 and dilepton production rate [19] in presence 
of dimension four electric and magnetic condensates were 
investigated some time ago. Also attempts were made to 
extract a nonperturbative electric screening mass from 
the gluon propagator [20 . In this work we calculate for 
the first time the structure of the screening masses in 
a very comprehensive way using OPE that provides us 
the nonperturbative information to the infrared sensitiv- 
ity of QCD, which in turn may help in constraining the 
various thermodynamic and spectral properties of high 
temperature QCD matter. 

In QCD the gluon polarization operator is not trans- 
verse in general, P^II^^y (P) 7^ 0. The most general ten- 
sorial structure of the in-medium gluon self-energy for 
an 0(3) invariant gauge fixing condition can be written 
as [21 

TT^^ (Cj,p) = TTiVl^ + nT^lu + ^mM^u + TT^^^, . (1) 
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We work here in covariant gauge and also omit the color 
indices for brevity. The Lorentz invariant single par- 
ticle energy and momentum are, respectively, given as 

uj = u ■ P and p = \J [u • P)^ — where is the four 
velocity of the heat bath and P = (po^p)- The projection 
operators are defined as [21] |22] 
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(color online) The nonperturbative ghost-gluon ver- 
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and Pf^ are transverse with respect to P^, while M^^ 



satisfies a weaker condition P^M^^ Py 
functions in ([T]) are extracted as 



T^fiu 1 TTt = - Pt TT^jy , 

Here tt^ and tt measure the deviation from transversal- 
ity. In high temperature perturbative QCD (pQCD), the 
transversality holds only in the temporal axial gauge and 
Feynman gauge [8 . In general the violation of transver- 
sality is however sub-leading in temperature in pQCD 
and one usually neglects TTm and tt. 

Now from ([T]) , the most general form of the gluon prop- 
agator = Po,/ii. (1 + "^fiu^o^fii^y^ follows as 



- p2_ 



2 [2 {p'-TTi) ir'p' 



[(^-lp2 _ ^) ^ + ^^M^, + (P2 - TT,) C^,] .(4) 



The Slavnov- Taylor identity (STI) in covariant gauge, 
P^V^yP"" = P^VQ^^yP"" - ^, leaves three independent 
components in Q. The most general form of the non- 
perturbative gluon propagator is 



= Pl,Di^Pl^Dt^M^yDm 



(5) 



in an obvious notation. 

The chromoelectric and chromomagnetic condensates 
are given by the second moment of the nonperturbative 
gluon propagator [17] 
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where Fa = — 1 and Nc is the number of color. Note 
that the frequency sum is only restricted to the lowest 
Matsubara mode (/cq = 0) in the spirit of the plane wave 
method [T^- This is equivalent to restricting oneself to 
the most dominant infrared singular sector. Also the 



lowest Matsubara mode excludes the explicit appearance 
of Dm in the gluonic condensates. Similarly, the ghost 
condensate is given as 



{f]Ur])^ = TFa 



(2^ 



G(0,/c) /c^ 



(7) 



The scalar where G is the nonperturbative ghost propagator. 
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FIG. 2. (color online) Gluon self-energy with gluon conden- 
sate (1^^ and 2^^ diagrams), ghost condensate (3^^ diagram) 
and ghost-gluon mixed condensates (4^^, 5^^, 6*^ plus their 
permutation) . 



Let us note some of the important points considered 
in our calculations: 

1. We work here to the effective order of ag in the 
sense that terms which are higher order in the cou- 
pling are related to terms of the order of as through 
the equations of motion. 

2. The gluon condensates are composite operators 
which do not correspond to any conserved currents 
and thus, are not renormalizable. Nonetheless, 
one can extract finite non-renormalizable contribu- 
tions by combining with other composite operators. 
The point is that under renormalization the glu- 
onc operators acquire admixtures of certain other 
operators, e.g., ghost. Therefore, we include the 
ghost-antighost condensate and two loop contribu- 
tions involving nonperturbative ghost-gluon vertex 
as shown in Fig. [l] This fixes uniquely the coeffi- 
cients of dimension four gluonic condensates f23] 
in the gluon self-energy as shown in Fig. [2] The 
moment of the nonperturbative ghost-gluon vertex 
in Fig. [T] can be obtained as. 



^3/ fabc^_ 



{2nr (27r) 



(8) 



3 



3. At zero temperature, the gluon propagator that in- 
cludes ah possible condensates up to mass dimen- 
sion four in Fig. [2] satisfies the STI where the 
ghost and mixed ghost-gluon condensates cancel 
the longitudinal terms generated by the gluon con- 
densate [15]. However, extending it to finite tem- 
perature we find that the STI is not obeyed [24] . 
This is not unexpected as it is also found in per- 
turbative [8] [21] as well as in nonperturbative cal- 
culations [20] . 

4. We note that the magnetic screening mass depends 
on the gauge fixing parameter ^. On the lattice 
one can measure various quantities by fixing the 
gauge. In the same spirit, we also intend, based on 
OPE with input from LQCD, to estimate screening 
masses for various gauge choices and compare with 
those lattice results. 

Now, the in-medium propagating modes can be written 
from (|4l) as 
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In OPE, the nonperturbative corrections to the polar- 
ization tensor are calculated by writing down the full 
Feynman diagrams and subtracting the equivalent per- 
turbative ones. The soft loop momenta are expanded in 
powers of external momenta and moments are identified 
with condensates as described above [25 . This is quite 
different from the Hard Thermal Loop approximation of 
pQCD where the polarization operator is saturated by 
the hard loop momenta (~ T). At finite temperature, 
the general expressions for scalar functions in the non- 
perturbative gluon self-energy are obtained by summing 
all the diagrams in Fig. [2] These functions are quite 
involved [24 and reveal a rich structure of thermal QCD 
in a nontrivial background. The screening masse^ can 
be extracted from the pole position of the propagator in 
the spacelike region po = 0, = — M^, where M is the 
relevant mass scale. The nonperturbative contribution 
to various scalar components of the gluon-self energy in 
the static limit po ^ can be obtained as 
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^ In Refs. ElIIO] a different prescription was used to extract screen- 
ing masses from gauge invariant correlators within a dimension- 
ally reduced QCD. This dimensional reduction works at a rather 
high temperature. A direct comparison of our non-perturbative 
prescription with those in Refs. [9. .10 may not be justfied. 
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«f(0,p)=0, 7r"''(0,p) 7^0 
where, 
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Here R = We = {2 + l) , Wb = ^ (38 + 9^ 

and Wm = (2 + 0- The value of a obtained here is 
the same as in Ref. [20]. Condensates appearing in the 
last two terms in tt^ are classical equations of motion for 
ghost and gluon fields so they vanish. The electric and 
magnetic screening masses are obtained from 
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To the perturbative order a^, yrp^^ (0,p) 
AiiasT'^, whereas Trf'^''^ (0,p) = 0. Solving ([l2| and ([l3|, 
we obtain the values of the screening masses as 
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For numerical evaluations of nonperturbative part in 
screening masses we use electric and magnetic conden- 
sates related to space (A^r) and timelike (A^-) plaquettes 
measured on lattice for pure SU (3) gauge theory [13 as 



TT 
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(15) 



where (S^)q is gluon condensate at T = and we take 
{Q'^)q/T^ = 2.5 and the critical temperature, = 260 
MeV, for pure SU{3) gauge theory. The perturbative 
piece is evaluated using the two loop running coupling 
constant, 



as (/i) 



47r 



2pi InL 



(16) 



with /3o = 11, /3i = 51, L = In (mVA^). We take Jl = 
27rT and A = 1.03Tc. 

We further note that there is no as dependence of the 



condensates in (15), i.e., they cannot be expanded as a 



power series with respect to a^, since they are based on 
non-perturtbative LQCD results [TJ- Due to the same 
reason, the parametric power counting hierarchy of scales 
mjj > rrim ~ g^T is obscure at low temperature where 
the condensates might provide plausible explanation for 
the magnetic mass. 
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FIG. 3. (color online) Temperature variation of electric 
(upper panel) and magnetic (lower panel) screening masses. 
The present investigation is represented by OPE with two 
different gauge fixing parameter ^ whereas LQCD data are 
represented by LAT-I [3^ and LAT-II 4 . PLO is perturbative 
leading order. 



The electric and magnetic screening masses so obtained 
are delineated in Fig.|3] We find that the nonperturbative 
contribution hmiD dominates over the perturbative lead- 
ing order (PLO) contribution for the temperature range 
we considered. The complete electric screening mass, in- 
cluding perturbative and nonperturbative contributions, 
falls short of lattice data but is still rather close to it. On 
the other hand, the magnetic screening is purely nonper- 
turbative in nature and agrees relatively well with lat- 
tice data. As seen the magnetic mass is dependent upon 
the gauge fixing parameter and we have chosen Landau 
(^ = 0) and Feynman (^ = 1) gauge. The weak gauge de- 
pendence found here is in agreement with that of Ref. [3] 
for a similar choice of gauge fixing. 

In summary, we have for the first time computed the 
nonperturbative contribution to both chromoelectric and 
chromomagnetic screening masses using OPE in a gluonic 
plasma. In particular, the magnetic screening mass is in 
relatively good agreement with the LQCD data whereas 
that of electric screening exhibits some discrepancy. The 
OPE electric screening mass is about 20% below the 
LQCD data points, which, however, show a rather large 
spread. The knowledge of these quantities sets the dy- 
namical length scale and provides us with the active de- 



grees of freedom in a hot QCD plasma. These results 
may be useful input to calculate various thermodynamic 
quantities, spectral properties and for the phenomenol- 
ogy of jet quenching [26 , quarkonium suppression ^7] in 
hot QCD matter produced in relativistic heavy- ion colli- 
sion experiments. 
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